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ABSTRACT 

We demonstrate the versatility of a dual imaging polarimeter working in tandem with a Lyot 
coronagraph and Adaptive Optics to suppress the highly static speckle noise pattern — the greatest 
hindrance to ground-based direct imaging of planets and disks around nearby stars. Using a double 
difference technique with the polarimetric data, we quantify the level of speckle suppression, and 
hence improved sensitivity, by placing an ensemble of artificial faint companions into real data, with 
given total brightness and polarization. For highly polarized sources within 0.5", we show that we 
achieve 3 to 4 magnitudes greater sensitivity through polarimetric speckle suppression than simply 
using a coronagraph coupled to a high-order Adaptive Optics system. Using such a polarimeter with 
a classical Lyot coronagraph at the 3.63-m AEOS telescope, we have obtained a 6.5a detection in 
the iJ-band of the 76 AU diameter circumstellar debris disk around the star HR 4796A. Our data 
represent the first definitive ground-based near-IR polarimetric image of the HR 4796A debris disk 
and clearly show the two outer ansae of the disk, evident in Hubble Space Telescope NICMOS/STIS 
imaging. Comparing our peak linearly polarized flux with the total intensity in the lobes as observed 
by NICMOS, we derive a lower limit to the fractional linear polarization of > 29% caused by dust 
grains in the disk. In addition, we fit simple morphological models of optically thin disks to our data 
allowing us to constrain the dust disk scale height (2.5^'g AU) and scattering asymmetry parameter 

(g=<cos#>= 0.20l'5o)- These values are consistent with several lines of evidence suggesting that the 
HR 4796A disk is dominated by a micron-sized dust population, and are indeed typical of disks in 
transition between those surrounding the Herbig Ae stars to those associated with Vega-like stars. 
Subject headings: instrumentation: adaptive optics — methods: data analysis — stars: individual 
(HR 4796A)techniques: image processing — 



1. INTRODUCTION 

Direct imaging surveys of the close-in (< 1000AU) 
environments of nearby stars for companions and 
circumstellar disks are comin g into maturit y and 
returning sp e ctacul a r results (lLafreniere et ail |2007|: 
Nielsen et all 120081: iKalas et all l2008t IMarois et al l 
20081 lOppenheimer etafl 12 008: Ma cintosh et alJ 120081 : 
Metchev fc Hillenbrandll2009D . Much of the technical ef- 
fort in high-contrast imaging centers around suppress- 
ing the overwhelming luminosity of the host star, and 
the residual speckle noise, largely caused by uncor- 
rected aberrations in the incoming wave front. This 
quasi-static source of noise, an especially bad prob- 
lem for ground-based efforts, is stable for minutes or 
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hours (|Hinklev et a l. 2007) and is the largest obsta- 
cle to direct detection of companions or circumstellar 
disks (iRacine et"aT1ll999l: ISivaramakrishnan et al.l 12002b 
IMarois et all 120031 : ISoummer et alJ l2007f ). Many au- 
thors have suggested useful techniques for direct sub- 
traction of speckle noise through image post processing 
(ISparks fe Fordll2002t IMarois et all 1200a iHinklev et all 
|2007| ). Another technique using a dual-channel imag- 
ing polarimeter can extract a polarized signal due 
to cirumstellar material (companion s or a disk) from 
the ar r ay of unpolarized s peckle s dKuhn et al.l l200ll: 
Potter! 120031 : iPerrin et all 12004 lOppenheimer et all 



2008) . The Lyot Project (lOppenheimer et al.ll2003l l200 



Sivaramakrishnan et al.l l2007f) employed this technique 
in addition to using a very high-order A daptive Optics 
(AO) system (|Roberts fc Nevmanll2002j) and a cla ssical 
Lyot coronagraph working at the diffraction limit (|LvotJ 
ll939tlMalbet1ll996l : ISivaramakrishnan et al.ll200ll ). 

In this work, we carry out a sensitivity analysis for 
speckle suppression through polarimetry, and quantify 
the level of suppression (Section 2). In Section 3 we 
demonstrate the power of this technique with a detection 
of the debris disk surrounding the nearby (d = 72.8 ± 1.7 
pc), young (8 ± 2 Myr) star HR 4796A (AO V, U=5.78 
mag,lJurJ(|1991h:lKoerner et~aTI (119981); lAugereau et all 



(1999); ISchneider et alJ (|1999f ): Ivan Leeuwenl (|2007h ). 
achieved in the ff-band from the ground. To our knowl- 
edge, this is the first high-contrast polarimetric near- 
IR image of the the HR 4796A disk obtained from a 
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ground-based observatory. This technique is especially 
well suited for regions of the image which are heav- 
ily dominated by speckle noise. Although circumstel- 
lar disk imaging has largely been performed using space- 
based observatori es (Krist et al. 2000; Stapc lfeldt et al.l 
l2004HKalas et al.ll2005t iGraham et al.ll2007l ). this paper 
demonstrates the power of speckle suppression and its 
ability to image circumstellar disks from ground-based 
observatories. 

2. SPECKLE SUPPRESSION THROUGH POLARIMETRIC 
OBSERVATIONS 

The advantages of differential imaging polarimetry for 
high-contrast observations of circumstellar disks have 
previo usly been discus s ed at length by several au- 
thors (iKuhri et al.l[200l[ iPotterl [2003 : iHaies et~aTll200fl 
iPerrin et al.ll2004L [2008f ). We briefly repeat here some 
basic concepts to establish consistency and notation. 
The polarization of light is usual ly represented by the 
Stoke s vector [I. O. U. U] (|Stoked Il852t IChandrasekharl 
Il946h . The usual astronomical convention is for the +Q 
direction to be oriented north-south, and +U northeast- 
southwest, with angles increasing counterclockwise from 
north to east. Linear polarization can also be expressed 
in terms of polarized intensity, (Q 2 + U 2 ) 1 / 2 , and posi- 
tion angle 6 = 1/2 arctan (U/Q). The normalized polar- 
ized intensity (Q 2 + U 2 + V 2 ) 1 ! 2 jl is referred to as the 
degree of polarization, polarization fraction, or percent 
polarization. Notation is not always consistent: some 
authors use P to refer to polarized intensity while others 
use it for degree of polarization. In this work, capital /, 
Q,U,V, and P will always refer to intensities (e.g., with 
units o f Janskies o r Jy a rcsec ~ 2 ), no t norm alized quan- 
tities. iTinbergenl (| 19961 ) and iKelleri (|2002ft provide ex- 
cellent introductions t o astronomical polarimetry, while 
lAdamson et al.l (|2005| ) summarize the recent state of the 
art. 

At its simplest level, a dual-channel differential imag- 
ing polarimeter consists of any device that splits an im- 
age into two orthogonal polarization states. This is fre- 
quently achieved through the use of a Wollaston prism, 
a two-element birefringent prism which separates an in- 
coming beam into two orthogonal polarization states, 
while at the same time introducing an angular deflec- 
tion between the two beams. The project described in 
this paper used such a configuration. A typical image 
showing the two displaced fields (left and right) result- 
ing from the beam deflection, and their directions of po- 
larization is shown in Figure [TJ A measurement of the 
Stokes Q parameter can be obtained through a difference 
of the left and right channels. Such a subtraction ("Dif- 
ference: +Q" in Figure[lJ largely eliminates the unpolar- 
ized speckle halo common to both left and right channels. 
However, to eliminate the bulk of the remaining aberra- 
tions (aberrations not common to both channels) that 
persist in this difference image, we obtain a second mea- 
surement by modulating the polarization states by 90°. 
Subtracting these in turn gives a —Q image. We mod- 
ulate the polarization through the use of Liquid Crystal 
Variable Retarders (LCVRs). This —Q image is obtained 
by swapping the positions of the polarization states and 
subtracting the two channels. After the subtraction, any 
astrophysical object will now possess negative counts in 
the image, but those non-common path aberrations will 



TABLE 1 

Table of HR 4796A Observations on UTC 2005 January 26 



Polarization Mode Filter Exp. Time No. of Images Total 

+Q H 120 s 3 360 s 
-Q H 120 s 3 360 s 
+U H 120 s 3 360 s 
Total 1080 s 

have the same sign and spatial characteristics present 
in the +Q image. Subtracting the — Q image from the 
+Q image (ideally) eliminates the non-common aberra- 
tions, leaving only the astrophysically interesting targets 
present. Different modulations of the polarization states 
can be used to obtain Stokes U and V. Since no pho- 
tons are lost in this process, the Stokes / image can be 
obtained by summing the left and right images of any 
polarization configuration. A schematic representation 
of the reduction process for a full polarimetric sequence 
is shown in Figure [TJ 

2.1. Observations 

Under very good observing conditions on 2005 Jan- 
uary 26 UTC at the 3.63 m Advanced Electro-Optical 
System (AEOS) telescope in Maui, we obtained three 
-ff-band (1.65 fim) coronagraphic polarimetric sequences 
(+Q, -Q, +U modes) of the star HR 4796A. The 
data were gathered using the Lyot Project coronagraph 
(iQppenheimer et al.ll2003L 120041 : ISivarama krishna n et al.1 
2007[) and the Kermit infrared camera ( Perrin et al.l 
20021) . The coronagrap h was a diffraction-limite d, clas- 
sical Lyot coronagraph (|Lvotlll939l : iMalbed 119961 ) with a 
455 mas diameter occulting mask and a hard-edged Lyot 
stop. The AEOS telescope is an altitude-azimuth de- 
sign, equipped wit h an AO system using a 9 41 actuator 
deformable mirror iRoberts fc Nevmanl (|2002fl . The total 
observing time for this dataset was 1080 s, comparable 
to the 1024 s for the lSchneider et alj (fl999l) F160W HST 
data. During the observations, the local Fried parame- 
ter, ro, a measure of the strength of turbulence in the 
atmosphere above the observatory, spanned the range of 
15 to 25 cm, indicating nearly ideal conditions for AO ob- 
servations at AEOS. Over the course of the observations, 
we obtained only +Q, -Q, and +U images because our 
retarders did not provide sufficient retardance to obtain 
a — U image. 

The polarimeter implemented in the Lyot Project for 
the data in this paper was unique in two regards. First, 
the Wollaston prism was located immediately after the 
Lyot pupil in the coronagraph. This post-pupil location 
is the correct location to minimize differential aberrations 
between the two beams. This setup is an improvement 
over other pol arimeters designed f or use in high-contrast 
imaging, e.g. IPerrin et al.l (|2008h . Second, the use of 
LCVRs as a polarization modulator is relatively rare for 
night-time polarimetry. The benefits of using LCVRs in- 
clude a great deal of increased flexibility in modulation, 
a lack of any image motions induced by rotating optics, 
and slightly faster modulation (although still not faster 
than the atmospheric timescales involved). Disadvan- 
tages of using retarders of this type include a potentially 
reduced wavefront quality. 

2.2. Polarimetric Dynamic Range 



An image of the HR4796A disk achieved through speckle suppression 
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Fig. 1. — A schematic illustration of the steps involved in the reduction of data taken with a dual-channel imaging polarimeter. A faint 
80% polarized artificial companion has been placed at the 3 o' clock position. The "LCVR Configuration 1" panel shows the two orthogonal 
polarization images produced onto the detector by the Wollaston prism, here labelled by the vertical and horizontal arrow symbols. A 
difference of the left and right images gives a measure of the Stokes Q amplitude and eliminates the strong speckle halo common to both 
the left and right images. Modulating the polarization states ("LCVR Configuration 2") switches the placement of the two orthogonal 
images, and the resulting subtraction of these images gives the — Q amplitude. A difference of the +Q and — Q images eliminates those 
aberrations not common to the left and right channels, while preserving any polarized signal (see text). The images have a square-root 
scaling, and the scale bars reflect the increasing sensitivity as the reduction progresses. 



Our goal is to quantify the gain in dynamic range 
achievable using the dual-imaging polarimetry technique. 
As a reference, we start by illustrating the dynamic range 
achieved on our Stokes I images without taking advan- 
tage of the polarimetric capabili ties. According to a 
technique we discus sed previously (jSoummer et al.ll2007t 
IHinklev et alll2007h . we have derived the magnitude dif- 
ference (dynamic range) between the occulted star and a 
5(7 point source as a function of position in the Stokes / 
images. The residual scattered light outside of our coro- 
nagraphic mask, usually in the form of highly persistent 
speckle noise, is the main limiting factor for detection of 
a point source (See Section 2.3.1). Local evaluation of 
the amplitude of this noise in turn determines the mini- 
mum brightness required for a 5ct point source detection. 
A radial curve of this sensitivity is shown in Figure [5] and 
labelled "Stokes-I Dynamic Range'.' 

All photometric values in this paper were calibrated to 
unocculted images of HR 4796A, directly prior to the oc- 
culted sequences. The raw data images were calibrated 
through standard dark current subtraction, application 
of bad pixel maps, and flat-fielding. The flat field im- 
ages were acquired using incandescent lamps each night. 



Also, binary star observations with well known orbits 
were observed to calibrate the pixel scale and image ro- 
tation fiducials. Prior to the subtraction, the images are 
rotated so that north is up in the image, east is to the 
left, and registered to each other using a cross correlation 
with subpixel accuracy. The data reduction techniq ue is 
discussed in gr eater detail in|Soummer et al.l (|2006t ) and 
Appendix A of lOppenheimer et al.l (|2008) 

To actually calculate the polarimetric sensitivity on the 
image, artificial point sources with varying total bright- 
ness and degree of polarization were placed radially into 
the mutually orthogonal left and right images in the orig- 
inal data. An example of this is shown in Figure [I] which 
shows an artificial 80% polarized (in the vertical direc- 
tion) faint point source with a given total brightness. 
The relative flux for the artificially placed companion in 
the left and right channels is given by: 

2^ = ^*2*21(1 -p cos 2(9) (1) 

F^ = ^f±(l+pco S 29), (2) 
where F to tai is the total flux of the object (F to tai — 
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Fig. 2. — Top panel: The three lines labelled "Stokes-7 Dynamic 
Range" show the faintest detectable H- band 5cr point source , rel- 
ative to the ho s t star , as described in Hinkley et al. (2007]) and 
ISoummer et alj 120071) . Each line indicates a co-addition of data 
with equivalent exposure time of 360 s, 720 s, and 1080 s. This 
sensitivity is limited by the bright, uncorrected speckle halo in the 
field-of-view. The collection of three lines labelled "100% Polari- 
metric Dynamic Range (AAfp [ ar )" show the faintest detectable 

100% polarized source using our dual-imaging polarimeter. The 
dual-imaging polarimeter removes a majority of the unpolarized 
speckle halo, and the polarimetric sensitivity increases with expo- 
sure time at a rate close to that expected for uncorrelated Gaus- 
sian noise. The lower panel shows the deviation of this improved 
polarimetric sensitivity, through co-adding two and three total po- 
larimetric sequences, from the expected ^/Atexp gain attributed 
to purely Gaussian noise. The inner grey region shows the size of 
our occulting mask, while the wider grey region shows where the 
suppression of the speckle noise was not effective. 

Ft + F<^>). These images were processed through the 
normal double differencing process schematically shown 
in FigureQ] and the total flux of these point sources in the 
original images were incrementally increased from zero 
until a 5cr detection was obtained in the double difference 
images. This threshold total intensity was recorded and is 
plotted in Figure [2] labelled "100% Polarimetric Dynamic 
Range (A¥ pok )." 

2.3. Results 

Figure [2] shows the speckle suppression achieved with 
the double differencing technique is most effective within 
1.75" and can achieve up to 4 magnitudes of improve- 
ment over the Stokes I within 0.5". This level of 
polarimetric suppr ession is comparable to other AO- 
based polarimeters ([Potter et alj2000HPe"rrin et alj2004 
2008), but results using this particular polarimeter on 
AEOS benefit from the ext remely high-order AO system 
(| Robe rts & Neyman 20 02T), and the optimized corona- 
graph ( Sivaram akrishn an et alj 12001) in the beam path 
before any of the polarimctry optics. 

Modelling point sources are especially useful in the con- 
text of the data discussed in the next section, since the 
HR 4796 A disk ansae are very near to point sources. 
Finding the required polarization for detection in our 



double-difference technique will thus help us to further 
constrain the value obtained directly (using published 
NICMOS Stokes I values) described in Section 3. These 
calculations using point sources can be directly applied 
to extended objects with resolved surface brightnesses. 
We have recasted the sensitivity results in terms of sur- 
face brightness (mJy/arcsec 2 ), and those values are listed 
on the right-hand axes of Figure [2j It should be noted, 
though, that in Figure [5] the brightness difference values 
(left-hand axis) indicate the total brightness of a point 
source, while the surface brightness (mJy/arcsec 2 ) val- 
ues reflect the brightness of the peak intensity of the 
point source. Nonetheless, our analysis for point sources 
translates over to extended sources since the key issue is 
the overall brightness of the source in comparison to the 
amplitude of the surrounding noise. 

2.3.1. Post Speckle Suppression: A Return to the Read 
Noise Limited Regime 

The uncorrected speckle noise in high contrast imag- 
ing data is due to highly static aberrations in the in- 
coming wave front that are not corrected by the AO 
system, or are induced "downstream" from the wave 
front sensor. Consequently, the resulting s peckle noise in 
the im ages is highly stable with time. As iHinklev et al.l 
(2007) demonstrate (esp. their Fig. 2), since the speckle 
noise pattern persists with sequential images, simple co- 
adding of data does not significantly improve the corona- 
graphic sensitivity, as would be expected for uncorrelated 
Gaussian-type noise. The speckles must be removed to 
gain improvements in sensitivity. iMarois et~aT1 ([2006h 
have performed speckle subtraction through image post- 
processing to greatly enhance their dynamic range, while 
IHinklev et alj ([20081 ) employ an instrument which uses 
the wavelength dependence of the speckles to disentan- 
gle them from a true astrophysical source. In this work, 
we use a dual-imaging polarimeter to subtract the unpo- 
larized speckle pattern from the images. 

Once the highly-static, highly time-correlated, speckle 
noise has been removed through polarimetry, the result- 
ing noise characteristics in the double difference image 
are distinctly similar to noise with Gaussian type prop- 
erties. This is reminiscent of the read noise dominated 
regime in which speckles are not the dominant source of 
noise. In this regime, normal Gaussian-like noise prop- 
erties of the image become applicable, and sensitivity to 
polarized sources should increase with the square root 
of exposure time. Indeed, our data show behavior quite 
close to this. The dashed and solid lines of the polari- 
metric dynamic range curves shown in Figure show 
the sensitivity with double and triple the effective ex- 
posure time of that represented by the dotted line. If 
the image noise that dictates the sensitivity is similar to 
Gaussian type noise, the expected gain in sensitivity is 
2.5 log V2 ~ 0.38 mag and 2.5 log V3 ~ 0.60 mag. In the 
lower panel, we show the deviation in the curves from 
this expected gain. These residual curves are within 0.1 
- 0.2 magnitudes of zero, consistent with a noise pattern 
with largely Gaussian properties. Moreover, the fact the 
the residual is consistently negative indicates that the 
speckle noise pattern is nearly, but not quite in the Gaus- 
sian regime. This is due to any residual speckle pattern 
that was not completely subtracted during the double 
difference process. 
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TABLE 2 

HR 4796A //-band Disk Morphology 



Study 




Ansae Separation 


PA (deg E of N) 


Schneider ct al. 


(1999) 


2.12 ±0.04" 


26.8 ±0.6° 


Schneider ct al. 


f2009) 


2.107 ±0.0045" 


27.01 ±0.16° 


This work 




2.10 ±0.05" 


27.5 ±2.5° 



3. DETECTION OF THE HR4796A DEBRIS DISK 

Using the method described in the previous section, 
we have obtained a modest, yet significant detection of 
the circumstellar disk around HR 4796A shown in Fig- 
ure [3] Although the full ringlike structure of the disk 
is not immediately evident in our image, we clearly de- 
tect the two extreme edges of the disk (ansae) at the 
6.5cr level above the residual image noise. We mea- 
sure a position angle of 27.5° ± 2.5°, in good agreement 
with t he 27°. 01 ± 0°.16 as measured bv [Schneide r et all 
(2009). The intensity in linearly polarized light (shown 
in figure [3]) was determined from our double-difference 
Q images (described schematically in Figure [1} as well 
as the trio of Stokes U images to construct a normal- 
ized Stokes P image: -P[ meai - = \/Q 2 + U 2 /I. We find a 
peak polarized flux density on the brighter (north) lobe 
of 7.4 mJy/arcsec 2 . Comparin g this to the pu b lished 
peak brightness (Stokes I) from Schneider et "all (1999) 
of 17 mJy/arcsec 2 , we derive a fractional polarization 
of 44 ± 5%, compara ble to values found in other debris 
disks, e.g. AU Mic ([Graham et al.ll2007l) . Expressing 
this as the 3<r lower limit, we state that the true polar- 
ization is greater than 29%. However, it should be noted 
these values are only lower limits to the true fractional 
linear polarization: recent Lyot Project data from the 
AEOS telescope suggest some Stokes V polarization in- 
duced by the AEOS telesc ope (|Qppenheimer et aT1l2008t 
lHarrington fc Kuhnl l2008) may be present, reducing the 
amount of linearly polarized flux in the Q and U modes. 
Since this particular observing sequence did not have the 
capabilities to measure Stokes V, we present our result 
merely as a lower limit. Given these issues, along with 
the relatively small telescope aperture, this is still a sig- 
nificant demonstration of the speckle removal technique. 

Shown in Figure [2] is the peak surface brightness 
(17 mJy/arcsec 2 ) of the bri ghter of the two disk ansae 
from I Schneider et alj ([1999) . Inspection of the plot sug- 
gests this corresponds to a polarization level of 70%, still 
consistent with our direct lower limit calculation of 44%. 
However, such a comparison may not be completely valid, 
as the dynamic range analysis assumed point sources for 
the sensitivity derivation while the disk ansae are more 
extended lobes. 

3.1. Morphological Models 

To constrain the nature of the debris disk around HR 
4796A, we model optically t hin disks assuming a Henyey- 
Greenstein phase function (|Henvev fc Greensteinl fl94ll) 
and a Raleigh-like variation of polarization with scatter- 
ing angle, which is suitable for disks with small grains 
or larger grain aggregrates (|Kimura et al.1 [20061 ). The 
model is a two-dimensional g eneralization of the o f the 
one-dimensional model used in I Graham et all (|2007[ ), ap- 
propriate for a solar sys t em zodiacal Henyey- Greenstein 
dust model (|Hondll985l : iGraham etld1l2007l ). We have 
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Fig. 3. — An //-band fractional polarization (Pjj, 



y Q 2 + U 2 /I) image for the star HR 4796A taken using the Lyot 
Project dual-imaging polarimeter at the AEOS telescope. The ex- 
treme lobes (ansae) of the circumstellar debris disk are detected at 
the 6.5ct level and show a lower limit fractional polarization of 44% 
(see text). The central black circle indicates the size of our coro- 
nagraphic mask, and the image has been smoothed with a 3 pixel 
(40 mas) smoothing function. The image i n the lower left corner i s 
an HST 1.6^m (F160W) image taken from [Schneider eTaLI llT99"9h . 



chosen to fit for the Henyey-Greenstein scattering asym- 
metry paramter g=<cos9> and the disk scale height. 
We fit for these two parameters, since both are intrinsi- 
cally related to the dust structure in the disk, and can 
most readily be constrained by the polarimetric data. A 
sample grid of models with scale heights ho = 12AU, 
25AU, and 50AU as well as Henyey-Greenstein param- 
eters 3=0.0, 0.3, and 0.6 is shown in Figure 2] to guide 
the eye. A value g = is completely isotropic scattering, 
while g — 0.6 signifies moderately strong forward scat- 
tering. A value of g = 0.3 is typical of Zodiacal dust and 
some debris disks. The model assumes a Gaussian ver- 
tical density distribution with an adjustable scale height 
based on the COBE model of zodiacal background light 
(|Kelsall et al.| [l998). Using cylindrical coordinates, we 
adopt the following density structure for the disk, 



n(r, z) = n 



J'o 



-0(z/h o ) 



(3) 



with a = -1 8 03, 8 = 4.973, 7 = 1.265 motivated 
by iKelsall et all (|1998l) . We have used the inner and 
outer radii (69 AU and 83 AU, respectively), inclination 
(14.12° from edge on) , and position angle (63.2°) from 
iSchneider et alj ([2009). Each thumbnail image in Fig- 
ure [5] was computed using the measured pixel scale of 
the Lyot Project's infrared camera (13.5 mas/pixel) and 
has been convolved with a PSF for the AEOS telescope 
at 1.65 /im. This PSF also reflects the reduced effective 
pupil diameter imposed by the Lyot mask in the Lyot 
Project coronagraph. During the fits, we also mask out 
the region covered by our coronagraphic mask. We gen- 
erated an ensemble of models varying ho between 2.5 and 
25AU and g between 0.0 and 0.7. 
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Fig. 4. — A matrix of fractional polarization images (same as Fig- 
ure [3j of our disk models for various morphologies. Each grid col- 
umn represents a different disk scale height (12, 25, or 50AU); each 
row shows a different value of the Henyey-Greenstein g=<cos8> 
parameter (0.0,0.3, and 0.6). The disk position angle, inc lination, 
inner/outer radii were taken from Schneider ct al. (2009). Gaus- 
sian noise has been added to match the data in Figure [3] 

3.2. model fit results 

The x 2 fitting procedure was performed on a pixel- 
by-pixel basis, normalized with the peak brightness in 
the model matched to the peak brightness in the polar- 
ization image. Ideally, such a fit should simultaneously 
use the polariz ed intensity and th e total intensity of 
the disk image dGraham et al.|[2007f ) in the fits to avoid 
degeneracies between dust properties a nd disk structure 
(jDuchene et al.ll2004t IPmte et al.ll2007| ). Given the rela- 
tively low signal-to-noise of this detection, and the lack 
of a detection in a total intensity image, we chose only 
to perform a \ 2 fit to the polarization image. Also, the 
data showing the fractional polarization may give a bet- 
ter handle on the disk scale height. In addition, the 
speckle noise pattern in highly corrected AO images can 
be well m odelled with a Rician Probability Distribution 
Funct ion (lAime fc Soummer|[200^ ; iFitzgerald &: Graham! 
120061: ISoummer et al.ll2007l) . so a model fit should take 
this fact into account. However, given that the speckle 
noise pattern has been significantly subtracted away, as 
discussed above, a \ 2 minimization employing a Gaus- 
sian noise variance is suitable. 

Our best \ 2 fit (x 2 — 1-1) is a model with dust 
disk scale height of 2. 51^'° AU and Henyey-Greenstein 
phase function of 0.20±;5q. These findings are consis- 
tent with several lines of evidence suggesting that HR 
4796A is dominated by a micron-sized dust population. 
Such a value of the asymmetry parameter is similar 
to models employing ISM-type dus t grain distributions 
(|Wood et al.ll2002t iKhn et al.lll994t) . with the dominant 
grain size lying between 0.1 and 1 /im. This population 
is similar th at of the (8±2 Myr Al V) 49 Cet disk as 
discussed bv lWahhai et alj (|2007[ ). with small (0.1 /im) 
grains between 30 and 60 AU from the star. On the other 
extrem e, in their analysis of HR 4796A, iDebes et al.l 
(2008) suggest that the minimum grain sizes for this disk 
are between 1 and 5 /im in size, with a dominant popula - 
tion of 1.4 ^m grains. Although the IDebes et al.l (2008) 



study employed the use of Tholins in their model fits, 
iKohier et all (|2008f ) find that the disk spectra can be 
well fit using less exotic compounds (amorphous silicates, 
amorphous carbonates, and water ice). In these models, 
1 (im grains give the best overall fit. Another powerful di- 
agnostic of the grain size distribution is the level of linear 
polarization measured in this work (44 ± 5%). Given the 
efficiency with which small grains can polarize starlight, 
our results are consistent with a micron-sized dust popu- 
lation leading to this relatively high level of polarization. 
The age and grain size di stribution for HR 47 96A sup- 
port the interpretation of iWahhai et alj (|2007t ). that it 
as well as the 49 Cet disk are representative of a class of 
disks that are "transitional" between the Herbig Ae class 
and Vega-like stars. 

Our best estimate and uncertainties on the disk scale 
height do not significantly constrain it below 12.5 AU at 
the 2cr level. Although several author s have measured 
gas s cale heights in disks such as f3 Pic (|Brandeker et al.l 
2004), to our knowledge this work is on e of the only 
measu rements of a dust disk scale height. iKenvon et al.l 
(119990 have simulated the structure of the HR 4796A 
disk, and predict a scale height between 0.3 and 0.6 AU, 
nearly consistent with our results. Our value for the 
phase function is similiar to the phase functi on values ob- 
taine d for other disks, in cluding Fomalhaut (IKalas et 
[2005]) and HD 141569A (|Clampin et al. I [20031 

4. CONCLUSIONS 

In this work, we have demonstrated and quanti- 
fied the speckle suppression achievable with a dual- 
channel imaging polarimeter working in conjunction with 
a high-order AO system and a classical Lyot corona- 
graph. Using a double subtraction technique in the 
data processing, we have demonstrated that a dual- 
channel imaging polarimeter can facilitate detection of 
polarized point sources within ~1" that are 10 — 100 
times fainter than the faintest (unpolarized) objects de- 
tectable by a coronagraph+AO system alone. These re- 
sults should be of inter est for future pola r imetry-based 
planet finding projects (iGisler et all 120041: iBeuzit et al.1 
120061: Macintosh etall 120061: iLucas et alj (2009f) . With 
this instrument, we report the first ground-based near-IR 
polarimetric detection of the disk surrounding the young 
star HR 4796 A Usin g Stokes I values from NICMOS 
(|Schneider et al.|[l999t ) to complement our data, we de- 
rive a lower limit to the fractional linear polarization of 
> 29% caused by dust grains in the disk. In addition, 
we fit simple morphological models of optically thin disks 
to our data allowing us to constrain the disk dust scale 
height (2.5jli*3 AU) and scattering asymmetry parame- 
ter (g=<cos9>= 0.20l'5o)- Given these values, and the 
relatively high level of polarization, these findings are 
consistent with several lines of evidence suggesting that 
the HR 4796A disk is dominated by a micron-sized dust 
population. 
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